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Abstract 

Limits on gravitational Cerenkov radiation by cosmic rays are obtained and used to constrain coefficients for Lorentz violation in 
the gravity sector associated with operators of even mass dimensions, including orientation-dependent effects. We use existing data 
from cosmic-ray telescopes to obtain conservative two-sided constraints on 80 distinct Lorentz-violating operators of dimensions 
four, six, and eight, along with conservative one-sided constraints on three others. Existing limits on the nine minimal operators 
at dimension four are improved by factors of up to a billion, while 74 of our explicit limits represent stringent hrst constraints on 
nonminimal operators. Prospects are discussed for future analyses incorporating effects of Lorentz violation in the matter sector, the 
role of gravitational Cerenkov radiation by high-energy photons, data from gravitational-wave observatories, the tired-light effect, 
and electromagnetic Cerenkov radiation by gravitons. 


7. Introduction. A century after its formulation. General Rela¬ 
tivity (GR) is established as a remarkably successful classical 
field theory of gravity. Extending GR into the quantum domain 
is widely believed to require modifications of one or more of its 
founding principles, and identifying experimental tests to con¬ 
firm this forms an interesting challenge. One central component 
of GR is local Lorentz invariance, which is symmetry under lo¬ 
cal rotations and boosts. Experimental tests of this invariance 
have undergone a renaissance in recent years [1], following the 
realization that tiny observable violations of Lorentz invariance 
could arise naturally in some quantum theories of gravity such 
as strings [2]. While impressive sensitivities to a broad range 
of possible violations in the matter sector have been achieved, 
searches for Lorentz violation in the gravity sector have been 
less extensive. In the present work, we obtain tight constraints 
on a large class of potential Lorentz-violating operators in the 
pure-gravity sector. 

Deviations from local Lorentz invariance in gravity can be 
studied using effective field theory, which offers a model- 
independent approach to describing Lorentz-violating effects 
arising in an underlying theory of quantum gravity [3]. Within 
this approach, the Lagrange density describing general Lorentz 
violation for pure gravity is a subset of the gravitational 
Standard-Model Extension (SME) consisting of the usual 
Einstein-Hilbert and cosmological-constant terms, along with 
a series of additional terms containing all possible Lorentz- 
violating operators. In each term, the Lorentz-violating opera¬ 
tor is formed from gravitational-field variables and is contracted 
with a coefficient controlling the magnitude of the effects. The 
Lorentz-violating operators can be organized in a series ac¬ 
cording to increasing mass dimension d in natural units, with 
the corresponding coefficients having mass dimensions A - d. 


Within the pure-gravity sector of this framework, most exper¬ 
imental studies [4—11] and theoretical investigations [12-17] 
have focused on minimal operators for Lorentz violation, which 
have mass dimension d - A. Some observational consequences 
of nonminimal operators of dimensions d - 5,6 are known 
[18], and experimental studies of nonrelativistic effects ofd-6 
operators on short-range gravity have recently been performed 
[19-21]. Eor reviews see, for example. Refs. [22-24]. Here, 
we obtain stringent conservative constraints on certain Lorentz- 
violating operators of even dimensions d > A, following from 
the observation of high-energy cosmic rays and the consequent 
limits on vacuum gravitational Cerenkov radiation. 

Electromagnetic Cerenkov radiation in ponderable media has 
been extensively studied since its discovery in the early 1930s 
[25, 26]. It arises when the velocity of a massive charged par¬ 
ticle exceeds the phase velocity of light in a medium, thereby 
rendering the particle unstable to radiation of Cerenkov light. 
In the presence of Lorentz violation, the vacuum acts like a 
refractive medium for particles with properties controlled by 
the coefficients for Lorentz violation [27]. Under these circum¬ 
stances, a particle travelling in a vacuum with velocity exceed¬ 
ing that of light can produce vacuum Cerenkov radiation, which 
continues until the particle loses enough energy to drop below 
light speed. The observation of high-energy particles of various 
species limits the existence of vacuum Cerenkov radiation and 
therefore places constraints on certain coefficients for Lorentz 
violation in the matter sector [28-40]. Any single coefficient 
constraint is normally one-sided because Cerenkov radiation is 
possible only for superluminal particles, which typically occurs 
for only one coefficient sign. 

Gravitational Cerenkov radiation is an analogous effect that 
is hypothesized to occur when the velocity of a particle exceeds 


Preprint submitted to Elsevier 


WHET 596, August 2015 




the phase velocity of gravity. In principle, this could occur in 
conventional GR in the presence of dust, gas, or other media, 
but the radiation rate is suppressed by two powers of the New¬ 
ton gravitational constant Gn and hence is negligible for prac¬ 
tical purposes [41^5]. However, in the presence of Lorentz 
violation, vacuum gravitational Cerenkov radiation suppressed 
by only one power of Gn can arise and would produce energy 
losses of particles travelling over astrophysical distances [46- 
49]. The observation of high-energy cosmic rays therefore con¬ 
strains certain coefficients for Lorentz violation in the gravity 
sector. In this work, we use observations of the energies and 
celestial positions of cosmic-ray events to obtain conservative 
limits on a large class of coefficients in the pure-gravity sec¬ 
tor, setting stringent first constraints on many coefficients and 
improving certain existing limits by factors of up to a billion. 

Observable effects on photon propagation arising from 
Lorentz violation involving operators of arbitrary d can be clas¬ 
sified in terms of anisotropy, dispersion, birefringence, and 
whether they affect vacuum propagation [50]. A comparable 
analysis of quadratic gravitational operators of arbitrary d re¬ 
veals that a similar classification holds also in the gravity sector 
[51]. We focus here on nonbirefringent vacuum effects involv¬ 
ing gravitational operators of arbitrary mass dimension, which 
can intuitively be viewed as certain components of a derivative- 
dependent effective metric s“^. We obtain wave solutions for 
this class of modifications to the Einstein field equations, de¬ 
rive expressions for the rates of vacuum gravitational Cerenkov 
radiation of scalars, fermions, and photons, and apply the re¬ 
sults to extract explicit conservative constraints on coefficients 
for Lorentz violation for even mass dimensions 4 < c/ < 8. 
Throughout this work, we use the notations and conventions of 
Ref. [3]. 

2. Lorentz-violating gravitational waves. The effective gravita¬ 
tional field theory containing Lorentz-violating operators of ar¬ 
bitrary mass dimensions [3] can be linearized to produce modi¬ 
fied Einstein equations relevant for weak-field gravity at leading 
order in coefficients for Lorentz violation [12, 18, 51]. Our fo¬ 
cus here is on perturbative modifications that can be written in 
the form 

Gfiv = S7tGn(Tm)hv + ( 1 ) 

where Gn is the Newton gravitational constant, {TM)fiv is the 
matter energy-momentum tensor, Rapys = ea/iKASySfiyR^^^''I is 
the double dual of the Riemann tensor, and G^y is the Einstein 
tensor. All expressions are understood to be linearized in a flat- 
spacetime background with Minkowski metric, g^y = rj^y + h^y. 
The operator s^y s Sy^ is formed as a sum of terms containing 
coefficients -for Lorentz violation and even powers 

of derivatives, 

^^y = (2) 

d 

with the sum understood to range over even values d > 4. The 
coefficients “ are constant and assumed to be small. 

The d - 4 coefficient s^y = (^^^)/jy appears in the minimal 


Lorentz-violating extension of GR and has been the subject of 
various experimental tests [4—11]. Nonrelativistic effects from 
some components of the second term (s'-^^)/jy^^ have recently 
been experimentally constrained as well [19-21]. 

The perturbative change to the field equations (1) preserves 
diffeomorphism symmetry even though the background coeffi¬ 
cients violate it. This can be understood as a con¬ 

sequence of the spontaneous breaking of diffeomorphism and 
Lorentz symmetry [52]. As a result, the usual counting of de¬ 
grees of freedom in the metric fluctuation /z^y holds, with four 
auxiliary components, four gauge components, and two phys¬ 
ical gravitational modes relevant to observable physics. The 
additional modes arising from the higher derivative powers ap¬ 
pear only at high energies that lie beyond the domain of va¬ 
lidity of effective field theory. Note that the form of Eq. (1) 
reveals that s^'' plays the role of a derivative-dependent shift of 
the metric, //^y —> rj^y - S/jy. The perturbation also includes 
Lorentz-invariant contributions, with effects governed by the 
trace pieces of the coefficients More general 

modifications to the Einstein equations can be countenanced 
and classified [3, 18,51], but exploring the implications of these 
lies outside our present scope. We remark in passing that odd 
powers of derivatives in the expression (2) are excluded by the 
requirements of hermiticity and diffeomorphism invariance. 

We seek solutions to the modified Einstein equations (1) rep¬ 
resenting perturbations of the conventional linearized gravita¬ 
tional waves propagating in the Minkowski vacuum. The wave 
solutions can readily be found at leading order in s^''. The con¬ 
ventional Einstein equations in vacuum 

Rf^v = 0 (3) 

hold at zeroth order, which implies that the modified Einstein 
equations at first order can be written in the form 

+T^)Ra^py ^ 0. (4) 

Note that working at first order in in this equation requires 
keeping both zeroth- and first-order contributions to the con¬ 
traction of the Riemann tensor with rf^, but only zeroth-order 
contributions to that with To fix the gauge, we adopt the 
modified Hilbert condition 

+ rp)hp, = 15 ^( 77 “^ + (5) 

and the traceless condition (rf^ +'s“l^)hap = 0 as natural choices, 
in light of the interpretation of the perturbation as a shift of the 
inverse metric. 

To find the wave solutions, it is convenient to convert to mo¬ 
mentum space via the ansatz 

(6) 

where 1“ is the 4-momentum of the gravitational wave and 
where as usual only the real part of the right-hand side is taken. 
This implies the replacement da —> Ha can be adopted in the 
definition (2) whenever l^y acts on h^y. Eor d - 4, the quadratic 
momentum dependence implies that the result tracks the con¬ 
ventional case modulo the deformation of the Minkowski met¬ 
ric. However, for d > 6 the corrections to the usual solutions 
involve higher powers of the 4-momentum P. 
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Using the ansatz (6), the modified Einstein equation (4) in 
the gauge (5) takes the form 

+^^)U/?A^v = 0. (7) 

The resulting dispersion relation for the gravitational waves, 

+ (8) 

suggests introducing an effective vacuum refractive index n - 

n(l) for gravitational waves, given by 

n2 = l-r^U^, (9) 

where la = la/lo- We emphasize that for d > 6 the refrac¬ 
tive index n depends on the momentum and hence the en¬ 
ergy of the gravitational wave and that for all d it receives 
direction-dependent Lorentz-violating contributions as well as 
both isotropic Lorentz-violating ones and Lorentz-invariant 
ones. The group velocity Vg can be obtained by differentiating 
l() with respect to I, yielding 

iu;i = 1 + 5 - 3)f-mna,... 

d 

( 10 ) 

At leading order in the wave vector la can be written in 
terms of the conventional wave vector in the form 

( 11 ) 

where - Sa^(l) is evaluated using Z®. The amplitude of the 
wave (6) can similarly be written in terms of the conventional 
plus and cross modes of GR. To obtain an explicit expression, 
it is convenient to work with trace-reversed quantities defined 
at first order by 

V = V + (12) 

where s rj^y - and = if'’ + In terms of the 
trace-reversed amplitude A^y, the gauge condition (5) then takes 
the form 

lair^-^+rP)Ap,^0. (13) 

This yields an expression for the trace-reversed amplitude A^y 
in terms of the conventional trace-reversed amplitude A®, 

(14) 

where the symmetrization on the indices ^,v includes a factor 
of 1/2. The result also shows that the graviton polarization ma¬ 
trix e^y appearing in the matrix elements for quantum processes 
with gravitons takes the form 

e,.=A(4«)-Vtv)X (15) 

where is the usual graviton polarization matrix and A is a 
normalization factor. 


3. Gravitational Cerenkov radiation. A particle of any species 
travelling faster than the phase velocity of gravity becomes un¬ 
stable to gravitational Cerenkov radiation. The differential rate 
to radiate a single graviton of momentum is given by 


czr = 


d^k 


dH 


\M\\2nf6\p-k-l), (16) 


2 -y/m2 H- jp- {27r)^2ko {2jif2lQ 


where is the incoming particle momentum obeying the dis¬ 
persion relation + ni^ for a particle of species w and 

mass m„, k^ — p^- Z^ is the outgoing particle momentum, and 
A1 is the matrix element for the decay in the quantum field the¬ 
ory. The integrated rate of energy loss is therefore given by 


dt 



d^k d^l 
(2n)^2ko (27r)3 


|M|2(2;r)V(p-k-Z). (17) 


Consider first the kinematical aspects of the rate of energy 
loss (17). Since the applications to follow involve m„ <K po, 
we can neglect contributions of order m„ times Lorentz viola¬ 
tion, which implies that in the delta function both Z and k are 
aligned with p at leading order. Performing the integrals over 
the outgoing momenta k and manipulating the remaining delta 
function yields 


dE 

dt 


8|^ 


mi, -I- jp- 


f 


dC 

(2;r)2|Zl 


|A1r(5(cos0- cos 0c). (18) 


where 0 is the angle between p and Z. The generalized Cerenkov 
angle 0c can be written as 


cos 0c = 


^ml+jp 1 


\P\ n(|Zt) 'AP\ \ [n(|zT)]2 


(19) 


where the convenient notation n(|Z|) s n(|Z|p) is used. The result 
(19) is a function of both \p\ and |Z|, and it encodes the Lorentz 
violation via the refractive index. 

The delta function in the integrand of the rate (18) governs 
physical properties of the gravitational Cerenkov radiation. Lor 
example, in the limit |Z| <K \p\ the integrand acquires contribu¬ 
tions only for radiation at the special Cerenkov angle 0 - 6c 
given by cos 0c = I In/3, where jS is the particle speed. This 
matches the well-known result for conventional electromag¬ 
netic Cerenkov radiation in a medium of refractive index n. The 
existence of a possible threshold velocity ySth below which no 
radiation occurs can also be seen. Lor example, in the above 
limit jSth = 1 /«. which again reproduces the classical result for 
photon Cerenkov radiation. The maximum angle of emission 
Oc, max = cos“'(l/n) occurs for an ultrarelativistic particle with 
y6 ^ 1 in this case. 

Lor calculational purposes, it is convenient to express the cor¬ 
rection to the refractive index arising from operators of dimen¬ 
sion d as 

n 1 - i (20) 

d 


where is a direction-dependent combination of coefficients 
for Lorentz violation given by 

P‘‘\l) = . . . W, (21) 
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and having mass dimension A - d. For the special case of an 
incoming photon or ultrarelativistic particle with —> 1, we can 
use the form (20) for n to obtain expressions for the Cerenkov 
angle 0c valid to leading order in Lorentz violation, 


cos 0c ~ l + 

d 

sin^Oc ^ ^(-1) 


d/2^{dhJld-4 


(d+2)/2^(dW^d-4 


1-^ 

V 1^/ 


1-^ 

\P\ 


( 22 ) 


The latter result reveals a cutoff at large graviton momentum I 
given by 


l^max ~ \p\' 


(23) 


The maximum momentum of a radiated graviton is thus approx¬ 
imately the momentum of the incoming particle, and this pro¬ 
vides an upper cutoff to the energy-loss integral (18). 

The explicit form of the integral (18) depends on the ma¬ 
trix element for the graviton emission, but dimensional anal¬ 
ysis shows that its basic structure is universal in the ultrarel¬ 
ativistic limit of interest here. At tree level, the emission of 
a graviton is proportional to the Newton gravitational constant 
Gff- Also, since the Cerenkov process is forbidden in conven¬ 
tional physics, the decay rate and the energy-loss rate must be 
controlled by the relevant dimensionless combination n - 1 of 
coefficients for Lorentz violation. The calculations of matrix 
elements performed below reveal that this dimensionless factor 
is (n - 1)^. Incorporating the various contributions to n from 
different d generates cross terms that complicate the calcula¬ 
tion, so for definiteness and simplicity we proceed here under 
the assumption that only a single value of d is of interest for a 
given analysis. The energy-loss rate then becomes proportional 
to The remaining dimensional factors must involve pow¬ 

ers of the incoming momentum \p\. The result of the integration 
(18) therefore takes the form 

^ = -F'^{d)GN{s^‘'^f\p\''‘'-\ (24) 

dt 

where F'"(d) is a dimensionless factor depending on d and on 
the flavor w of the particle emitting the gravitational Cerenkov 
radiation. The time of flight of the particle from its initial en¬ 
ergy Ei to a final energy Ef is then 


t = 


T'^'id) 

GM(iW)2 


rp2d~5 


j2d-5 


(25) 


where = {2d - 5)!F'^{d) is another dimensionless factor. 

To obtain explicit expressions for F'^{d) and T'^id), we must 
consider the matrix elements for specific processes. Here, we 
discuss in turn the cases where the incoming particle is a scalar, 
a photon, and a fermion. 


where e - -y/j^ is the vierbein determinant. The Feynman rule 
for the scalar-scalar-graviton vertex is i VTbTrGivC^v, where 

= -p^lkv + J0^v(p“ka + ml). (27) 

The second term has no effect in practice as it generates a trace 
of the graviton polarization in the matrix element and therefore 
vanishes for physical states. The squared matrix element for the 
tree-level process takes the form 

|M|2 = 167rG^C^,C^^e^’'e“^, (28) 

where with r - +,x are the two physical graviton polariza¬ 
tion modes contained in the matrix (15). Summing over these 
modes and inserting the Cerenkov angle, we find 

|M|2 = l6nGN{s^‘^^f\f^-^ {\pf - 2\l\\p\^ + \P\^W) . (29) 

Upon integration, the dimension-dependent factor T'^id) ap¬ 
pearing in Eq. (25) is found for massive scalars w s 0 to be 

T‘<’{d)^\{d-2){d-3). (30) 

In the special limit with only operators of mass dimension d - A 
and only isotropic effects, the results match those obtained in 
Ref. [46]. 

Photons radiating gravitons. Next, consider gravitational 
Cerenkov radiation from a high-energy photon. The electro¬ 
magnetic part of the Einstein-Maxwell Lagrange density is 

X = (31) 

where F^^ - df^Ay - dyA^ is the usual electromagnetic field 
strength. In TT gauge, the vierbein determinant contributes 
only at nonlinear order in the metric fluctuation and can be ne¬ 
glected here. Using the standard plane-wave ansatz for the pho¬ 
ton, the Eeynman rule for the photon-photon-graviton vertex is 
found to be i VTbTiGivC^vip, where 

G^yXp — 2tlpv{pdkp -f- k^Pp) "t" 2Pv{t]Xpkp -t rjppkd) 

+ jkp{ri2yPp + l]pyPA) - {pakaitlApTIpv + tlAvdpp)- (32) 

The square of the matrix element for the tree-level diagram is 
then 

\Mf - 2nGNCpyApCapys^,e“e'^efd/ef, (33) 

where the physical photon polarization vectors are ^ with 
s - 1,2 and repeated indices r, s, and t indicate sums over 
polarizations. An extra factor of 1/2 has been incorporated as 
usual for the sum over incoming photon polarizations. 

In evaluating the result (33), the condition p^k'’CpyAp = 0 
effectively implements the replacement —> rf^. It is con¬ 

venient to choose the 3 axis along I and the 2 axis along the 
component of p perpendicular to I, which gives 

Pi -0, P 2 - -PQ sin 0, Ppf - -po^o(l - n cos 0), (34) 


Scalars radiating gravitons. We first consider gravitational 
Cerenkov radiation from a hypothetical real massive scalar min¬ 
imally coupled to gravity via the Lagrange density 

£ = -\eg^'’dp(j)dy^ - jemlcj)^, (26) 


and leads to the identities 

PpPvi^^Y“i^'^ya = PiiPviY'Y^ie^y^ = \p\^ sin^ 0, 
(e")pv(e")^^ = (e%y(e^Y^ - 2, 
PpPviY)^'' = -\pf sin^ 0, PpPyie^)^'’ = 0. (35) 


4 



Using these results simplifies the squared matrix element to the 
form 


the coefficients for Lorentz violation must satisfy the one-sided 
constraint 


|M|2 = 

X (if? - A\i^\p\ + + 2\$\p\^ - \p\^) . 


(36) 


S^‘^\p) = ■ ■ -Pa,^. 




(44) 


Integration reveals that the dimension-dependent factor T^id) 
in Eq. (25) is 


THd) = 


{d- \){d-2){d-3){2d-3) 
4(4t/4 - 2%(P + 65d^ - 62d + 27) 


(37) 


for photons w = y radiating gravitons. Note that a key dif¬ 
ference between the massive-particle and photon cases is that 
photons are always above threshold for gravitational Cerenkov 
radiation. 


Fermions radiating gravitons. With the notation and conven¬ 
tions used in Eq. (12) of Ref. [3], the Lagrange density describ¬ 
ing the minimal gravitational coupling of a relativistic fermion 
can be written as 


£ = jiee^a^y^Df.i// + h.c., (38) 

where the gravitational degrees of freedom appear in the vier- 
bein and in the covariant derivative via the spin connec¬ 
tion. The leading contribution to the fermion-fermion-graviton 
vertex can be obtained by expanding the vierbein and noting 
that the spin connection contributes only at higher order. The 
Eeynman rule for the vertex takes the form i where 

Cty^\y,{Py + K). (39) 

The squared matrix element then becomes 

|M|2 = 87rG^ ^ (fi(p)C;f,M(k)e^’')(M(k)C^^M(p)e“^) , (40) 

spins 

where a factor of 1/2 is included for the usual average over the 
incoming fermion spins. Using the methodology developed for 
the photon case along with the standard projection 

^u{p)®u{p) = y p (41) 

spins 


for a relativistic fermion, we find 

|M|2 = 2,nGN{s^‘^^f\p\\t^‘‘-^ 

X (-lip + 3|;^|ip - 4|;^pi| + 2|;^^) . (42) 

The dimension-dependent factor T'^id) in Eq. (25) resulting 
from the integration for w = ij/K 


^ ^ (d-2)(d-3)(2d-3) 

^ A{2d^-ld + 9) 


(43) 


Note that this differs from the scalar result (30) due to an addi¬ 
tional term arising from the sum over spins. 


4. Constraints. According to the above results, the observation 
of a cosmic ray of species w arriving at the Earth with energy Ef 
after travelling a distance L along the direction p^ implies that 


This reveals that high-energy particles originating at large dis¬ 
tances offer the best bounds. 

The species dependence in the bound (44) is encoded en¬ 
tirely in the factor T'^id), which is given for scalars, photons, 
and fermions in Eqs. (30), (37), and (43), respectively. Note 
that these factors are finite for any finite d, and they vanish only 
for physically irrelevant values d < A. Note also that they imply 
enhancements in radiated power for increasing particle spin. As 
a result, for fixed values of £/ and L, photons yield more sen¬ 
sitive bounds than fermions, and fermions more sensitive ones 
than scalars. The improved sensitivity increases with increas¬ 
ing d. Eor example, the ratio R = T^(d) : T'^id) : T‘^{d) is 
R ^ 0.6 : 0.8 : 1 for d - A, but changes to R ^ 0.05 : 0.5 : 1 for 
d - b and becomes R —» d^^ \ d^^ 12 \ 1 at large d. 

To gain some initial intuition about the implications of the 
bound (44), consider the conservative scenario of a heavy nu¬ 
cleus travelling a distance L ^ 10 Mpc ^ 10^^ GeV ' from 
a nearby active galactic nucleus and impacting the Earth with 
an observed cosmic-ray energy of about 100 EeV. Assuming 
the gravitational Cerenkov radiation occurs from a partonic 
fermion in the nucleus carrying about 10* GeV of the total 
cosmic-ray energy and taking the factor 'F'^(d) to be of order 
1 for simplicity, we find constraints on combinations of coef¬ 
ficients for Lorentz violation of dimension A — d of approxi¬ 
mate order GeV"^^'^. Although only a crude estimate, 

this serves to reveal the quality of constraints from gravita¬ 
tional Cerenkov radiation. Eor example, bounds on some d - A 
Lorentz-violating operators should exceed by several orders of 
magnitude the various existing sensitivities, which are of order 
10“^-10“^° on dimensionless coefficients in the gravity sector 
[4-12]. Similarly, limits for the case d - 6 should reach 10“^* 
GeV“^ or so, representing stringent first constraints on this class 
of nonminimal coefficients in the gravity sector. 

Repeating the above crude estimate but replacing the imping¬ 
ing cosmic ray with a high-energy photon reveals that gravi¬ 
tational Cerenkov radiation from photons generically provides 
weaker constraints. For example, even the observation of an 
ultra-high-energy gamma ray at 100 TeV would imply a sen¬ 
sitivity to Lorentz-violating operators of dimension d reduced 
by a factor of compared to cosmic rays. This factor 

overwhelms any possible gain in sensitivity from greater pho¬ 
ton propagation distances L, even for the most favorable case 
with d - A and for photons originating at cosmological dis¬ 
tances. We therefore focus here on constraints from gravita¬ 
tional Cerenkov radiation by cosmic rays, deferring further con¬ 
sideration of high-energy photons to the discussion section be¬ 
low. 

To obtain more definitive constraints, information about the 
direction of travel of the cosmic rays is required, in addition to 
their energy and distance of travel. Since cosmic rays impinge 
upon the Earth from many directions on the celestial sphere. 
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Observatory 

events 

Emax (EeV) 

Ref 

AGASA 

22 

213 

[57, 58] 

Ely’s Eye 

1 

320 

[59] 

Haverah Park 

13 

159 

[57, 60] 

HiRes 

11 

127 

[61] 

Pierre Auger 

136 

127 

[62] 

SUGAR 

31 

197 

[57, 63] 

Telescope Array 

60 

162 

[64] 

Volcano Ranch 

2 

139 

[57, 65] 

Yakutsk 

23 

160 

[57, 66] 


Table 1: Cosmic-ray events and maximum energies used in this work. 


d j 

Lower bound 

Coefficient 

Upper bound 

4 0 

-3 X 10-14 < 

■'on 


4 1 

-1 X 10-1^ < 
-8 X 10-14 < 
-7 X 10-14 < 

■'lO 

Re i^41 
Im^|l 

<7x 10-14 

< 8 X 10-14 

<9x 10-14 

4 2 

-7 X 10-14 < 
-7 X 10-14 < 
-5 X 10-14 < 

-6 X 10-14 < 
-7 X 10-14 < 

“^20 

Re 41 
Im^jl 

Re ^22 

Im '^22 

< 1 X lO-iii 
<7x 10-14 
<8x 10-14 

< 8 X 10-14 

<7x 10-14 


Table 2: Conservative constraints on dimensionless coefficients 


d j 


Lower bound 

Coefficient 

Upper 

bound 


^(6) 

“^00 

< 

2 

X 

10 

-31 

-6 

X 

10 

-30 

< 

^U) 

< 

1 

X 

10 

-29 

-6 

X 

10 

-30 

< 

Re 

< 

7 

X 

10 

-30 

-8 

X 

10- 

-30 

< 

Im 5 !®! 

< 

5 

X 

10 

-30 

-1 

X 

10- 

-29 

< 


< 

1 

X 

10 

-29 

-7 

X 

10 

-30 

< 

Re 4? 

< 

7 

X 

10 

-30 

-9 

X 

10 

-30 

< 

Im^jl 

< 

6 

X 

10 

-30 

-9 

X 

10- 

-30 

< 

Re^l* 

< 

6 

X 

10 

-30 

-8 

X 

10- 

-30 

< 

Im 

< 

6 

X 

10 

-30 

-1 

X 

10 

-29 

< 

-(6) 

< 

8 

X 

10 

-30 

-8 

X 

10 

-30 

< 

Re ^3®! 

< 

7 

X 

10 

-30 

-6 

X 

10- 

-30 

< 
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it is natural to work with coefficients for Lorentz violation ex¬ 
pressed in spherical coordinates rather than cartesian ones [50]. 
The combination s*"'' of coefficients, which appears in the re¬ 
fractive index (20) and is given in terms of coefficients for 
Lorentz violation by Eq. (21), is an observer scalar and hence 
can be expanded in terms of spherical harmonics as 

,W(p) = ^ (45) 

jm 

where jm are the usual angular quantum numbers, subject here 
to the restriction that j is even and j < d - 2. At each fixed d, 
there are (d - 1)^ independent spherical coefficients. The result 
(44) then becomes a constraint on the spherical coefficients 

To obtain conservative bounds, we suppose that the cosmic 
ray primary is a nucleus of atomic weight N. Most high-energy 
cosmic rays are believed to be protons {N - 1), but some may 
be light nuclei or even heavy nuclei such as iron {N - 56) [53]. 
We also assume that any gravitational radiation is emitted by 
one of the fermionic partons in the nucleus, which carries a 
fraction r of the cosmic-ray energy observed at the Earth. A 
conservative estimate is r = 10% [46]. We therefore take the 
energy £/ as £/ = rEq^jN ^ £e/560, with the factor T'^(d) 
in the bound (44) identified as 'F'^(d) in Eq. (43). This con¬ 
servative estimate therefore represents a reduction by a factor 
(^/^)"'-5/2 ^ ^^Qd-5/2 j-jjg effective energy and hence of the 
bounds. The acceleration sites of cosmic rays are believed to 
be extragalactic, including possibly supermassive black holes 
in active galactic nuclei [54]. The nearest of these lies at a dis¬ 
tance of a few Mpc, which offers a sense of the minimum value 
of the distance L. The distance is limited by spallation of the 
cosmic ray on photons in the cosmic microwave background 
[55]. Approximately 50% of protons and iron nuclei are be¬ 
lieved to survive at distances of 100 Mpc, while for lighter nu¬ 
clei the analogous distance is 20 Mpc [56]. Eor definiteness, we 
take L ^ 10 Mpc. 

Numerous collaborations have published data for the ener¬ 
gies and angular positions of observed cosmic rays. Since 
higher-energy events provide greater sensitivity to coefficients 
for Lorentz violation, we restrict attention here to events with 
energies above 60 EeV. Table 1 provides some information 
about 299 observed events of this type. The first column lists 
the cosmic-ray observatory, the second shows the number of 
published events above 60 EeV used in this analysis, the third 
gives the maximum observed energy £max, and the final col¬ 
umn provides the reference. To obtain numerical constraints, 
we adopt the modified simplex method of linear programming 
[67] detailed in Ref. [68] in the context of bounds on Lorentz 
violation from nongravitational Cerenkov radiation by neutri¬ 
nos. In the present instance, the available dataset of 299 events 
is sufficiently large to place constraints on all coefficients for 
fixed d - A, d - b, and c/ = 8 in turn. In principle, higher 
values of d could also be considered, and additional cosmic-ray 
data for energies below 60 EeV could be included as well. 

Although the bound (44) is one sided for each cosmic ray, 
the dependence on the direction of travel and the plethora of 
data across much of the celestial sphere mean that independent 


6 



J 

8 


8 


8 


8 


8 


j 

Lower bound 

Coefficient 

Upper 

bound 

0 

-7 

X 

10 


< 

-(X) 

“^00 


1 

-1 

X 

10 

-45 

< 

^10 

< 

1 

X 

10 

-45 


-9 

X 

10 

-46 

< 

Re 

< 

8 

X 

10 

-46 


-9 

X 

10 

-46 

< 


< 

9 

X 

10 

-46 

2 

-9 

X 

10 

-46 

< 


< 

1 

X 

10 

-45 


-1 

X 

10 

-45 

< 

Re if/ 

< 

8 

X 

10 

-46 


-8 

X 

10- 

-46 

< 

Im 

< 

9 

X 

10 

-46 


-1 

X 

10- 

-45 

< 

Re If/ 

< 

9 

X 

10 

-46 


-1 

X 

10- 

-45 

< 

T -(s) 

Im S 22 

< 

9 

X 

10 

-46 

3 

-1 

X 

10- 

-45 

< 

-(«) 

“^30 

< 

1 

X 

10 

-45 


-1 

X 

10- 

-45 

< 

Re 

< 

8 

X 

10 

-46 


-9 

X 

10 - 

-46 

< 

Im^j' 

< 

9 

X 

10 

-46 


-8 

X 

10 - 

-46 

< 

Re 4 ) 

< 

9 

X 

10 

-46 


-9 

X 

10 

-46 

< 

Im^32' 

< 

8 

X 

10 

-46 


-8 

X 

10 

-46 

< 

Re ^3^ 

< 

1 

X 

10 

-45 


-1 

X 

10 

-45 

< 

Im^33> 

< 

1 

X 

10 

-45 

4 

-1 

X 

10 

-45 

< 

*40 

< 

1 

X 

10 

-45 


-6 

X 

10 

-46 

< 


< 

1 

X 

10 

-45 


-8 

X 

10 

-46 

< 

Im 

< 

1 

X 

10 - 

-45 


-8 

X 

10 

-46 

< 

Re^z 

< 

1 

X 

10 

-45 


-6 

X 

10 - 

-46 

< 

Im 

< 

1 

X 

10 

-45 


-7 

X 

10 - 

-46 

< 


< 

1 

X 

10 

-45 


-8 

X 

10 - 

-46 

< 

Im * 

< 

8 

X 

10 

-46 


-1 

X 

10 

-45 

< 

Re If/ 

< 

8 

X 

10 

-46 


-9 

X 

10 - 

-46 

< 

Im 

< 

6 

X 

10 

-46 

5 

-1 

X 

10 - 

-45 

< 

*50 „ 

< 

1 

X 

10 

-45 


-8 

X 

10 - 

-46 

< 

Re 7^5*^ 

< 

1 

X 

10 

-45 


-8 

X 

10 - 

-46 

< 

Imlf/ 

< 

7 

X 

10 

-46 


-9 

X 

10 

-46 

< 

Re 4 

< 

9 

X 

10 

-46 


-8 

X 

10 

-46 

< 

Im7<| 

< 

8 

X 

10 

-46 


-1 

X 

10 

-45 

< 

Re 4 

< 

7 

X 

10 

-46 


-6 

X 

10 

-46 

< 

Imif/ 

< 

1 

X 

10 

-45 


-9 

X 

10 

-46 

< 

Re 4 ) 

< 

1 

X 

10 

-45 


-8 

X 

10 - 

-46 

< 

Im 

< 

8 

X 

10 

-46 


-8 

X 

10 - 

-46 

< 

Re 4 

< 

1 

X 

10 

-45 


-8 

X 

10 - 

-46 

< 

Im 7 <! 

< 

1 

X 

10 

-45 

6 

-1 

X 

10 - 

-45 

< 


< 

2 

X 

10 

-45 


-8 

X 

10 - 

-46 

< 

Re 4 

< 

1 

X 

10 

-45 


-7 

X 

10 - 

-46 

< 

Im # 

< 

9 

X 

10 

-46 


-1 

X 

10 - 

-45 

< 

r 4 ' 

< 

6 

X 

10 

-46 


-6 

X 

10 

-46 

< 

ImJ 

< 

1 

X 

10 

-45 


-7 

X 

10 

-46 

< 

Re 

< 

1 

X 

10 

-45 


-7 

X 

10 

-46 

< 


< 

8 

X 

10 

-46 


-8 

X 

10 

-46 

< 

Re 

< 

1 

X 

10 

-45 


-9 

X 

10 

-46 

< 

Im #) 

< 

8 

X 

10 

-46 


-8 

X 

10 - 

-46 

< 

Re#) 

“^65 

< 

9 

X 

10 

-46 


-8 

X 

10 - 

-46 

< 


< 

9 

X 

10 

-46 


-1 

X 

10 - 

-45 

< 


< 

9 

X 

10 

-46 


-7 

X 

10 - 

-46 

< 

1 . 4 ' 

< 

1 

X 

10 

-45 


Table 4: Conservative constraints on coefficients in GeV 

jm 


two-sided constraints are implied for almost all spherical coef¬ 
ficients at each fixed d. The exception is the isotropic coef¬ 
ficient which produces orientation-independent effects and 
hence can only be constrained on one side. For definiteness, we 
perform two independent analyses at each d = 4,6,8. One as¬ 
sumes only the isotropic coefficient is nonzero and yields 
a single one-sided constraint. The second assumes a purely 
anisotropic model allowing all the coefficients with j ^ 0 
to be simultaneously nonzero and yields d(d - 2) independent 
two-sided constraints. Note that the spherical coefficients are 
complex when m + Q,so their real and imaginary parts must be 
treated as independent for this analysis. 

Tables 2, 3, and 4 contain the resulting constraints on the 
spherical coefficients s®, and respectively, obtained 
using cosmic-ray data and reported in the Sun-centered frame 
[69]. The initial two columns provide the values of d and j, 
while the corresponding spherical coefficients are listed in the 
third column. The numerical lower and upper bounds are given 
in the second and fourth columns, respectively, with units as 
specified in the table captions. For coefficients with d - 4, the 
results in Table 2 represent improvements of factors of a thou¬ 
sand to a billion over existing maximal sensitivities obtained via 
direct laboratory measurements [1]. Note that the connection 
between the spherical coefficients and the usual cartesian 
ones is given by equations analogous to Eq. (130) of Ref. [40]. 
For coefficients with c/ = 6,8, the results in Tables 3 and 4 are 
the first constraints in the literature. For d - 6, they comple¬ 
ment the constraints on other independent coefficients obtained 
in experiments testing short-range gravity [19-21]. The reader 
is reminded that the constraints in tables 2, 3, and 4 are con¬ 
servative: if the cosmic rays are assumed to be protons instead 
of iron nuclei and if the full proton energy is assumed available 
for gravitational Cerenkov radiation, then the displayed bounds 
for d = 4, 5, and 6 would be sharpened by additional factors of 
more than 10"^, 10^, and 10'^, respectively, even for the same 
propagation distance L. 


5. Discussion. In this work, we have derived properties of grav¬ 
itational waves in the presence of a class of Lorentz-violating 
operators of arbitrary d, used the results to derive energy losses 
from gravitational Cerenkov radiation, and performed an analy¬ 
sis of existing cosmic-ray observations to extract constraints on 
a variety of coefficients for Lorentz violation with d = 4,6,8. 
With the exception of the bound on the isotropic coefficient 
all the measurements reported here are the first constraints 
obtained from gravitational Cerenkov radiation on the corre¬ 
sponding Lorentz-violating terms, and none of the coefficients 
for t/ = 6 or 8 have previously been constrained in the literature. 

The constraints in Tables 2, 3, and 4 are obtained using cos¬ 
mic rays rather than photons because the former are observed at 
much higher energies than the latter. Nonetheless, the absence 
of gravitational Cerenkov radiation from high-energy photons 
does in principle contain additional information. Indeed, in a 
more general treatment, each species w would provide distinct 
constraints because the particle itself experiences Lorentz vio¬ 
lation that is flavor dependent [27]. 


7 



To illustrate this, consider the Lorentz-violating vacuum dis¬ 
persion relation 

- P^-ml ^0, (46) 


where n„ = n„(p) is a refractive index for the particle of species 
w, with niH, = 0 if the particle is a photon. The motion of the 
particle then follows a geodesic in a pseudo-Finsler spacetime 
[70, 71]. An example would be a refractive index for a massive 
fermion given by [40, 72] 

nl^l+lZ^p^PH, (47) 

in analogy with the refractive index (9) for gravity, where (\v 
is a momentum-dependent coefficient having expansion of the 
form (2). For d - A this reduces to a special case of the matter 
sector of the minimal SME with a single fermion flavor [27], 
with Zv - Zv being a constant dimensionless coefficient for 
Lorentz violation. A similar result holds for neutrinos [39] and 
for photons [50], where the corresponding constant dimension¬ 
less coefficient is conventionally denoted by Repeating 

the analysis of gravitational Cerenkov radiation generates an 
energy loss given by Eq. (18) as before, but with the vacuum 
Cerenkov angle 9c given instead by 


cos 6c 


^Jml + Z[nU\p\-\i\)? 


\P\ 



ml+p^ 


im 


nw(\p\)n(\l\) 

Mih - \i\)f 
[«(lf)]2 . 

\ [nMp\-\i\)Z 


(48) 


where as before the arguments of the refractive index are under¬ 
stood to be oriented along p. This expression reveals that the 
presence and rate of gravitational Cerenkov radiation depends 
on ratios of the refractive indices for the particle and the gravi¬ 
tational waves. Depending on the relative magnitudes of the co¬ 
efficients for Lorentz violation for the particle and the graviton, 
gravitational Cerenkov radiation may occur in a given direction 
of travel for only one sign of the correction to the graviton 
refractive index, for either sign, or not at all. Note that the final 
term depends only on the particle refractive index, contributing 
only when depends on momentum and hence only for d > A, 
implying the particle is experiencing nonminimal Lorentz vio¬ 
lation. Note also that the above expression reduces to the result 
(19) in the limit n„, —> 1, as expected. 

Inserting the gravitational Cerenkov angle (48) into the in¬ 
tegral (18) for the energy loss and limiting attention to fixed 
d must by dimensional arguments produce a result of the gen¬ 
eral form (24) but now involving a linear combination of the 
quadratic terms (s*'^^)^, and (1 -n„,)^, where 1 -n„ is 

given in terms of matter coefficients for Lorentz violation by an 
expression analogous to that for in Eq. (21). It follows that 
a detailed analysis of observations of high-energy particles, in¬ 
cluding photons and neutrinos, would yield constraints on dis¬ 
tinct combinations of coefficients for Lorentz violation. How¬ 
ever, photons yield weaker constraints than those from cosmic 


rays because observed cosmic-ray energies exceed the highest 
photon energies by about a millionfold, which becomes scaled 
by the power c/-5/2 in extracting a bound. A similar conclusion 
holds for neutrinos. Note that any observed but unexplained 
absence of ultra-high-energy particles such as neutrinos or pho¬ 
tons could in principle be attributed to Lorentz-violating vac¬ 
uum Cerenkov radiation, including the gravitational Cerenkov 
radiation considered here. A complete analysis along these var¬ 
ious lines would be of interest but lies beyond our present scope. 

The expression (48) for 9c also reveals that if the Lorentz vi¬ 
olation is minimal, so that all Lorentz-violating operators have 
mass dimension d - A and both n and n„ are momentum in¬ 
dependent, then the properties of the gravitational Cerenkov ra¬ 
diation are determined by differences of SME coefficients for 
Lorentz violation in the gravity and matter sectors. This leads 
to an interesting relation between distinct measurements, as fol¬ 
lows. In this scenario, the rate of energy loss for a massive 
fermion undergoing gravitational Cerenkov radiation is gov¬ 
erned by the combination + 2Zf, while that for a radiating 
photon is governed by Moreover, due to the free¬ 

dom to choose coordinates without changing the physics, all 
nongravitational searches for Lorentz violation involving these 
species must involve the combination 2c[f^ + (kfZ^ [73]. As a 
consequence, if analyses yield measurements Mi of + 2^, 
M 2 of - {kcY^, and M 3 of TZZ -i- (kf)“^, then the three 
measurements must satisfy the condition 


Ml- M 2 - Mj.^ 0. (49) 

A relation of this type, relevant for searches for CPT violation 
with neutral-meson oscillations, has previously been inferred 
among SME coefficients in the quark sector [74]. 

The components of the coefficients (T('^))t'''“i -“<(-4 constrained 
in this work can also be measured in other ways. Lor exam¬ 
ple, all the corresponding Lorentz-violating operators for d> A 
are dispersive and so can in principle be measured using pulse- 
spread data from a terrestrial gravitational-wave observatory 
such as the Advanced Laser Interferometer Gravitational-Wave 
Observatory (LIGO) [75] and Advanced Virgo [76] or from a 
space observatory such as the proposed Evolved Laser Inter¬ 
ferometer Space Antenna (eLISA) [77], assuming gravitational 
waves are indeed detected. The sensitivity of dispersion mea¬ 
surements of depends on the ratio of the observed pulse 
width to the source distance L, which typically leads to weaker 
constraints than those from gravitational Cerenkov radiation. 
However, dispersion limits are distinct and unique in detail be¬ 
cause they involve Lorentz violation in the electron and pho¬ 
ton sectors by virtue of the laser interactions with the mirrors. 
Moreover, dispersive measurements involve the group veloc¬ 
ity and hence are sensitive to the gradient of the refractive in¬ 
dex instead of the refractive index itself, so the correspond¬ 
ing measurements constrain different combinations of Lorentz- 
violating operators when more than one value of d is incorpo¬ 
rated. 

Another interesting open issue is the prospects for inde¬ 
pendent two-sided bounds on the isotropic coefficients 
which cannot be obtained via gravitational Cerenkov radiation. 



Gravitational-wave observatories such as LIGO and eLISA are 
uniquely suited to place dispersion limits on these coefficients 
for operators of dimension d > A. However, no dispersion oc¬ 
curs for the minimal isotropic coefficient which makes its 
two-sided measurement challenging. One possibility in prin¬ 
ciple would be to compare the time of flight of gravitational 
waves to that of light or neutrinos emitted from the same source. 
This has the disadvantage of requiring simultaneous observa¬ 
tion using different techniques. Other options already used to 
place constraints on include analyses of data from Grav¬ 
ity Probe B [9] and from pulsar timing [10]. Methods such as 
the study of orbital decay rates of binary systems [78] have the 
potential to provide interesting sensitivities to as well. 

The derivations in this work may also have implications for 
other ideas. For example, the presence of Lorentz violation in 
quantum electrodynamics can induce triple photon splitting in 
the vacuum and hence reduce the frequency of light as a func¬ 
tion of distance travelled, suggesting the possibility of modifi¬ 
cations to the usual interpretation of the observed cosmological 
redshift [79]. A similar possibility is implied by the expres¬ 
sions (24), (25), and (37) for the energy loss of a photon due to 
gravitational Cerenkov radiation. These ideas are conceptually 
akin to ‘tired-light’ models [80], which are strongly constrained 
by the direct observation of time dilation associated with cos¬ 
mological redshift [81]. However, the energy losses (24) here 
have distinctive frequency dependence, and in principle they 
might only be perturbative or only affect part of the observed 
redshifts, perhaps such as those associated with supernova stud¬ 
ies of dark energy. A complete discussion of these possibilities 
would require analysis of gravitational-wave propagation in a 
cosmological background instead of the static Minkowski back¬ 
ground adopted here. Nonetheless, the analysis in the present 
work suffices to provide simple estimates of the possible scale 
of the effects, as follows. 

For an astrophysical source at small redshift z defined in 
terms of the photon energies as usual by z -i- 1 =£’;/£/, the 
luminosity distance Ll can be written as Ll ~ (z + 
where //q - 1.5 x 10 GeV is the Hubble constant [82]. Di¬ 
rectly expressing the time of flight (25) as a distance L in terms 
of z gives 


L 


GNis^Ef^-^ 


((z+l)2‘'-5 


1 )- 


(50) 


Comparing Ll and L reveals that the potential contribution of 
gravitational Cerenkov radiation to the observed cosmological 
redshift is primarily governed by the dimensionless ratio 


R = 




10‘*£f(n - l)^ 


(51) 


where £, is measured in GeV and values ^ > 1 represent sub¬ 
stantial effects. This suggests that redshift modifications from 
gravitational Cerenkov radiation are negligible for most prac¬ 
tical purposes. For example, for the optical frequencies ^100- 
900 nm typically studied in the spectra of type-la supernovae, 
the energy factor is of order 10“^^^ GeV^ or smaller, so a 
value /? > 1 would require n - 1 > 10'°, which is well out¬ 
side the perturbative regime. Moreover, the energy dependence 


implies that the effect varies by orders of magnitude over an 
observed spectrum, which for large values of n - 1 would dis¬ 
tort spectra beyond observed limits. The requirement of per¬ 
turbative n - 1 evidently restricts substantial redshift effects to 
high-energy photons. However, it remains conceivable that a 
detailed analysis along the above lines could extract additional 
constraints on from precision cosmological measurements. 

In the above, we consider gravitational Cerenkov radiation 
by photons. However, Lorentz violation in the pure-gravity 
sector can also cause electromagnetic Cerenkov radiation by 
gravitons, corresponding to graviton decay. Discussion of this 
process is lacking in the literature. The form of the Einstein- 
Maxwell Lagrange density (31) reveals that at leading order 
this process involves two-photon emission, being governed by 
a photon-photon-graviton-graviton vertex. The corresponding 
amplitude is proportional to Gn, and hence on dimensional 
grounds the power loss of the graviton takes the form 

^ = -F'(d)Gl(n - 1)VT, (52) 

where F'(d) is a dimensionless factor depending on d arising 
from integration of the matrix element, k is the power of the 
dimensionless combination n - 1 emerging from the matrix ele¬ 
ment, and \l\ is the magnitude of the graviton momentum. This 
result implies that the frequency spectrum of gravity waves de¬ 
tected by a gravitational-wave observatory on or near the Earth 
is distorted and downshifted by Lorentz violation. However, the 
effect is far below observational levels in practice, both because 
the power loss (52) is proportional to and because the en¬ 
ergy of typical gravitational waves is expected to be tiny. Eor 
example, a gravitational wave in the LIGO band at frequency 
100 Hz originating in our galaxy experiences a negligible fre¬ 
quency shift 6v ss 10“^^°(n - I)*’ Hz. Eor similar reasons, gravi¬ 
ton Cerenkov decay into other particle species is negligible as 
well. 

The results in this work complement those obtained in tests 
of short-range gravity [19-21] and thereby improve the cover¬ 
age of sensitivities to coefficients for Lorentz violation in the 
gravity sector. Exploring the remaining coefficients for even d 
and the coefficients for odd d, all of which are birefringent, is 
an interesting open problem for future research. 

This work was supported in part by the United States Depart¬ 
ment of Energy under grant number DE-SC0010120 and by the 
Indiana University Center for Spacetime Symmetries. 
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